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We present recently defined jet-observables for hadron-hadron dijet pro- 
| duction, which are designed to reconcile the seemingly conflicting theo- 

. retical requirement of globalness, which makes it possible to resum them 

(automatically) at NLL accuracy and the limited experimental reach of 
detectors, so that they are measurable at the Tevatron and at the LHC. 
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' 1. Introduction 

Qh, 

ft, ■ Event shapes and jet-rates are infrared and collinear (IRC) safe observ- 

QJ \ ables, which describe the energy and momentum flow of the final state. They 

constitute an ideal compromise between simplicity and sensitivity to prop- 
erties of QCD radiation. They provide then a wealth of information, e.g. 
in measurements of the coupling a s and its renormalisation group running, 
in cross checks/measurements of the values of the colour factors of QCD 
and, most importantly, in studies of the connection between parton-level 
(the perturbative (PT) description of quarks and gluons) and hadron-level 
(the real), for a review see 

IRC-safety ensures that event-shape distributions can be computed within 
perturbation theory, however in the more exclusive phase space region where 
perturbative radiation is suppressed (conventionally associated to almost 
vanishing values of the observable, V <C 1) large logarithmic corrections 
need to be resummed to all orders. More specifically, given an event-shape 
V, a function of all secondary final state momenta {ki, . . . , k n } and of the 
set of recoiling Born momenta {p}, the probability of "constrained events" 
i. e. V({p}, k\ . . . k n ) < v has a divergent PT expansion for v — > 

£(«) = Pvoh{V < v) = 1 + ^ Rn, m a s n ln m v + . . . , 

m<2n 
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i. e. there is a soft and collinear divergence [~» \nv] for each emitted gluon. 
Today's state-of-the art accuracy accounts for all Leading (LL) and Next- 
to-Leading Logarithms (NLL) as follows: 

S(«) = exp{Lgi(a s L) + g 2 (a s L) + . . .} , L = In - . (1) 

LL NLL 

Furthermore resummations are matched to fixed order results at NLO. 



2. Basics of resummation 

The main ingredient in resummations is factorisation. Usually two steps 
are needed to obtain a resummed prediction. The first task is to exploit an- 
gular ordering to reduce n-parton matrix elements to a QED-like factorised 
form, i. e. to the product of independent emissions from the hard scattering 
partons only. Schematically, e.g. for e + e~ — ► 2 jets this can be written as 

, . 1 -|-r , . 1 -pT- a s C 'f dE d9 

w pP (k u ...,k n ) = ^[[w pP (k i ) ~ — ■[[—— — — , (2) 

i=l ' i=l 

with corrections which contribute beyond NLL to £(«)• The second step 
requires some analytical understanding of the observable's behaviour in the 
presence of soft-collinear emissions in order to factorise its definition via 
Mellin transforms, i. g. for the thrust T in e + e~ — >2 jets one has (Q = ^) 

^ 1=1 1=1 
By combining eq. (J2J) and eq. © one obtains the resummed answer 



f dv VT [ f d9 dE a s (E6)C F ( 

Lit) / e exp / — — e 2C ? — 1 

v 1 J 2mu JOE vr \ 



(4) 



where the effects due to the running of the coupling have been included. 

Despite the simplicity of the ideas underlying resummations, some tech- 
nicalities cannot be avoided. Indeed, at NLL accuracy, care is needed to 
treat properly the emission of hard collinear or soft large-angle gluons and 
the inclusive gluon splittings. Furthermore, it is known that for some observ- 
ables, when seeking NLL accuracy, the multi-parton matrix-element cannot 
be factorised as in eq. ©, notably this is the case for non-global observ- 
ables [2*. In a similar way some observables, such as jet rates in Jade-like 
algorithms 0, do not exponentiate, so that it is not possible to write the 
distribution in the form in eq. (^Q). Also, the observable's factorisation itself 
can be non trivial especially for multi-jet observables like the 3-jet limit of 
thrust minor, which involves 5 integral transforms .41. 
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3. Automated resummed predictions 



The observation that the origin of logarithmic enhancement in jet-shape 
distributions is always the same (i. e. it is due to the radiation of soft- 
collinear gluons) and that the approximations which lead to resummed ex- 
pressions such as eq. (@J are very similar for different observables, it was 
natural to investigate the possibility of developing a general and automated 
approach to resummation, which does not require analytical treatment of 
observable specific multiple emission effects via Mellin/Fourier transforms. 

Such an approach is feasible, and has been implemented in the com- 
puter code CAESAR (Computer Automated Expert Semi- Analytical Resum- 
mer). jS] This program differs from usual Monte Carlos whose task is gener- 
ally the numerical evaluation of integrals via generation of random events. 
CAESAR instead is based on a generic master resummation formula, pre- 
sented and derived in detail in [5]. This formula applies to a well-defined 
class of observables, whose requirements can be summarised as follows. 

• Given a Born event, when just one soft emission k is radiated collinear 
to the Born parton p£, the observable should behave as 



where, for each hard leg I, ag,bg,dg are some numbers and ge(4>) is a 
regular function parameterising the azimuthal dependence. 

• The observable should be continuously global, meaning that it should 
be sensitive to emissions everywhere and the transverse momentum 
dependence should be uniform (a\ = . . . = a n = a). 

• The observable should be recursive infrared and collinear (r-IRC) safe, 
meaning that the addition of emissions which are much softer or more 
collinear should not drastically change the value of the observable. This 
new concept was introduced and illustrated in detail [3]. 

Given the general master formula, together with it's applicability conditions, 
CAESAR works as an expert system, which in a first step establishes whether 
the observable is within its scope. It's next task is to determine the inputs 
needed for the evaluation of the master formula, basically the coefficients 
and functions ai,bi,di, ge(4>) together with a NLL correction function T 
which accounts for effects due to multiple emissions jSJ E] • Finally as a last 
trivial step, CAESAR evaluates the master formula (integrating over different 
Born configurations when necessary). 

Notice that the first two steps are critical: they require high precision 
arithmetic [7j to take asymptotic (soft & collinear) limits and they follow 
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methods of "Experimental Mathematics" [H] to validate or falsify hypothe- 
sis. 

Currently implemented processes are 2 & 3 jets in e + e~, [1+1] & [1+2] 
jets in DIS, Drell-Yan + 1 jet, hadron-hadron dijet events. Results as ob- 
tained with CAESAR have been tested against all known NLL global event 
shape resummations in e + e~ , DIS and DY.IN hadronic dijet events no event- 
shape resummations had ever been carried out, so there was no standard 
definitions of event shapes, as in e + e~ and DIS. In the following we will 
therefore show how to define observables at hadronic colliders which recon- 
cile the seemingly conflicting theoretical requirement of globalness and the 
limited experimental reach of detectors. Indeed in hadron-hadron collisions 
the presence of radiation from the initial state, together with the limited 
coverage close to the beam region (usually expressed in terms of a maximum 
rapidity 77 max ) makes this a critical issue. 



4. Observables in hadron-hadron dijet events 

4-1. Directly global observables 

Measurements for this class of observables include ideally all particles in 
the event, i.e. the maximum rapidity ry max is taken as large as experimentally 
possible. One then defines variants of the usual e + e~ observables , e. g. 

• The transverse thrust 

Hi W±i ■ n T \ 



Tt = max ■ 



where q±i are the transverse momenta with respect to the Beam axis 
and ut defines the transverse thrust axis. 

The thrust minor 

El n ou t I 



here q° ut are the momentum-components out of the event-plane (i. e. 
that containing the beam axis and the n^-axis) • 

NLL resummations are then valid as long as hx(l/v) < (a + 6 m in)^max 0- 
Explicit results for the differential distributions show that most of the data 
lie in this region |10j . 



4-2. Observables with recoil term 

The drawback of observables presented in Sec. 14.11 is that they require 
measurements quite far in the forward region. While this is possible at the 
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Tevatron (where r/ max ~ 3.5), only calorimeter information is available in 
forward regions, whose resolution is quite low. A way to circumvent this is to 
define measurements which explicitly select a restricted, central region and 
to modify the observable's definition so as to make it indirectly sensitive to 
emissions in the unobserved phase space region. This sensitivity is typically 
achieved exploiting recoil effects. The idea is essentially that momentum 
conservation ensures that the vectorial sum of transverse momenta in the 
central region, exactly balances the vectorial sum of transverse momenta in 
the unobserved forward regions. The definition of this type of observables 
proceeds then as follows. One first selects a central region C (e. g. taking all 
particles with r\ < r/Q ~ 1), one then defines a central, non-global observable, 
e. g. a central thrust minor 



^Ei^i, q±, c = E l&<l. (6) 

and adds to it a recoil term, e. g. 



'-m,C — ^ 



1 



iec 



(7) 



Q±,e 

The recoil enhanced thrust minor is then given by 

T m ,n = T m fi + TZ±,c • (8) 

These observables are often named indirectly global observables, since they 
are indirectly sensitive of emissions in the unobserved region through recoil 
effects. The drawback of these observables is that in the region of extremely 
small V the NLL function J 7 , parameterising multiple emission effects UJ, 
has a divergence. The presence of this divergence is well understood, and is 
not specific of indirectly global observables. A divergence occurs whenever, 
due to cancellation between contributions from different emissions, an ob- 
servable can be very small in the presence of radiation. The origin of the 
divergence can be understood with simple arguments: the master resum- 
mation formula assumes that the mechanism responsible for keeping the 
value of the observable small is a LL Sudakov effect. However, in the phase 
space region of very small V, it is more likely that V is small because of 
bi-dimensional cancellations in TZ±£- The NLL function T cannot compen- 
sate a wrong LL behaviour and manifests this with a divergence. Despite 
the presence of the divergence, explicit results show that most of the dis- 
tribution is in a region which is under control of NLL resummations jlUj . 
What makes these observables challenging experimentally, is the accuracy 
with which the recoil term can be measured, indeed 7Z±£ is subject to big 
cancellations between almost equal and opposite large transverse momenta. 
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4-3. Observables with exponentially suppressed forward terms 

A variant of directly global observables is to define a class of observables 
with exponentially suppressed forward terms. 

One introduces the mean transverse-energy weighted rapidity rjc and 
total transverse momentum of a given central region C 

VC = "Fj E Vi Q±i ' Q±,c = ^2<1-Li, (9) 

^ L ' c iec i&c 

and with particle in the forward region one defines an exponentially sup- 
pressed forward term 

^—T^-E^i^ 1 ^ 1 , (10) 

then e. g. the thrust minor with exponential forward suppression is 

T m ,£ = T m fi + £c • (11) 

These observables are intended to have some of the better features of both 
purely global and indirectly global observables. Theoretical predictions are 
not affected by divergences and from an experimental side these observables 
are optimal in that there is no need for a fine resolution in rapidity and az- 
imuth in the forward region, and calorimeter information should be enough 
to determine the forward term Eq. 



5. Final considerations 

Despite the fact that observables with exponentially suppressed forward 
terms seem to be more suitable in many respects, we point out that observ- 
ables in different classes have complementary sensitivities to perturbative 
and non-perturbative (NP) radiation, so that the simultaneous study of a 
number of observables within each class is a powerful tool to investigate 
properties of QCD radiation, e. g. for 

• studies of underlying event, since the forward sensitivity (to beam- 
fragmentation) can be specifically tuned, suppressing it for purely PT 
studies, or deliberately enhancing it when studies NP effects. This al- 
lows one to test quantitative, as well as qualitative features of existing 
models, e. g. the sensitivity of the underlying event to the partonic 
channel (i.e. qq -> qq, qg qg or gg -> gg). 

• studies of hadronization corrections in multi-jet events, in particular 
for tests of power-corrections beyond the "Feynman Tube model" 
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• studies of the non-trivial quantum evolution of colour, indeed the novel 
perturbative QCD colour evolution structure that arises in events with 
4-jet topology ^5] has never been investigated experimentally before. 

Given the much more challenging environment in hadronic collisions, com- 
pared to e + e~ and DIS, we believe that only an automated approach can 
make such a programme feasible. A number of other possible measurements 
are given in ^U], while results for the automated analysis and resummation 
of a large number of observables are available from |http : / / qcd-caesar . org | 
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